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ABSTRACT 

We confirm our recent prediction of the “pitchfork” foreground signature in power spectra of high- 
redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all 
baselines. This is due to the inherent response of a wide-field instrument and is characterized by 
enhanced power from foreground emission in Fourier modes adjacent to those considered to be the 
most sensitive to the cosmological H i signal. In our recent paper, many signatures from the simulation 
that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this 
key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity 
through the coherent averaging of 12 independent snapshots with identical instrument settings and 
provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field 
effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple 
models for antenna apertures that have been proposed for future instruments such as the Hydrogen 
Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground 
leakage from the pitchfork by ^ 40 dB relative to the MWA and significantly increase the likelihood of 
cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum. 

Subject headings: cosmology: observations — dark ages, reionization, first stars — large-scale structure 
of universe — methods: statistical — radio continuum: galaxies — techniques: 
interferometric 
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1. INTRODUCTION 

The epoch of reionization (EoR) commenced follow¬ 
ing the formation of the first stars and galaxies. It is 
characterized by a period of non-linear growth of mat¬ 
ter density perturbations and astrophysical evolution in 
the universe’s history. The detection of redshifted 21 cm 
radiation of H i from this epoch is one of the most promis¬ 
ing probes of t he evolution of large-scale structure dur¬ 
ing this epoch (ISunvaev fc Zeldovi cbllI972HScott fc R^ 


200a) . 


I99rt iMadau et al.l Il997t iTozzi et al.l l200flt llliev et al 


Sensitive instruments such as the Square Kilometre 
Array (SKA), which will be capable of providing di¬ 
rect imaging of redshifted H i, are yet to become op¬ 
erational. In the meantime, the Hydrogen Epoch of 
Reionization Arraj0 (HERA), currently under devel¬ 
opment, will be much more advanced in its capability 
to detect and place definitive constraints on the reion¬ 
ization epoch relative to current instr uments such as 
the Murchison Widefield Ar r ay (MWA; [ Lonsdal e et al.l 
120091 : iBowman et al.l 1201,11 iTingav e t al.l 120 1, H). the 
Low Erequency Array (LOFAR.: Ivan Haarlem et al.l 
EqU), and the Precision Arr ay for Probing th e 
Epoch of Reionization (PAPER; iParsons et al.l I2OIOD . 
which have only enough sens itivity for a statis- 
tical detection of the signal (|Bowman et al.l 120061 : 


http://reionization.org/ 
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Parsons et alJl2012at iBeardsley et ^120131: iDillon et alJ 

2013t iThvagaraian et al.ll2013HPober et al.l[2014ll. 

The primary challenge to the detection of cosmological 
H I from the EoR arises from continuum emission from 
Galactic and extragalactic foreground objects, which 
is ^ 10"^ stronger than the desired signal. However, 
the inherent differences in spatial isotropy and spectral 
smoothness can be exploited to extract the cosmolog- 
ical signal from f oregro u nd co nta mination (see, e.g., 
Di Matte o et al.l 120021 12004 iFurlanetto fc I3ri^ 

2004 [Morales fc Hewitt] l2004 iZaldarriaga et af 


iSantos et al.l 120?^ Furlanetto et al. _ 20061 


McQuinn et al.~l200d iMorales et al.l 120061 iWang et alJ 


200i IGleser et aLI 12008^ Thus 


a detailed char- 
has become 


essential (lAli et al. 

120081: IBernarc 

i et al. 

120091 

20101; IBowman et al. 

12009: iLiu et al.l 

20091 

2014ayH: 

Datta et al.l 120101: Liu & TegmarkI 20111: IGhosh et al.l 

20121 Morales et al.l 120121: iParsons 

et al.l 

12012bl 

Trott et al.l 20121: 

Dillon et al.l 2014 Dillon et al.l 

2014 iPober et al.ll20l4 IThvagaraian et 

alJI201.4l201,'ih. 


referred to as Paper I) used instrument and foreground 
models, for the first time with full sky coverage, in or¬ 
der to simulate actual FoR experiments more accurately 
than previous studies. Surprisingly, we found that fore¬ 
ground emission outside the primary beam field of view 
caused the most significant contamination of the Fourier 
modes considered the most sensitive for detecting the 
cosmological H i signal in delay spectrum based analy¬ 
ses. This is the result of the interplay between foreground 
emission, particularly diffuse Galactic emission, and the 
wide-held properties typical of FoR instruments. Our 
simulations predicted that the delay spectra from the 
MWA and other experiments should exhibit a character¬ 
istic “pitchfork” appearance, with local maxima near the 
horizon delay limits, in addition to at the primary lobe 
region. 

Careful antenna aperture design can signihcantly mit¬ 
igate this contamination. Optimal weighting of contami¬ 
nated Fourier modes may be required to extract the sig¬ 
nal with maximum sensitivity. Thus, detailed knowledge 
of foreground signatures is key for the design and analysis 
choices of future instruments such as HFRA and SKA. 

In Paper I, we verihed the general features of our sim¬ 
ulations against MWA observations, but were unable to 
conhrm the pitchfork prediction due to an insufficient 
sensitivity in the small amount of data analyzed. Here, 
we use deeper MWA data to conhrm with a high sig- 
nihcance the presence of key pitchfork characteristics 
of wide-held measurements predicted in the preceding 
study. 

Section [2] is an overview of the role of wide-held mea¬ 
surements in the delay spectral domain and the predicted 
pitchfork signature. Section [3] describes the analysis of 
MWA data used to improve the dynamic range of the de¬ 
lay spectra. Section |4] describes the results and conhrms 
the presence of the predicted wide-held effects. Section [5] 
underscores their impact on aperture design. Section |6] 
summarizes our hndings. 

2. WIDE-FIELD EFFECTS IN DELAY SPECTRUM 

Paper I describes in detail the effects of wide-held mea¬ 
surements as seen in the delay spectra of interferometer 


visibilities. We list the relevant equations and give a brief 
overview of the wide-held signatures predicted therein. 

The delay spectru m for a baseline ve c tor, b, is 
(jParsons et al.l l2012allbl : IThvagaraian et aP 1201,4 Pa¬ 
per I): 

Vi,iT) = J Vtif)Wif)e^^^f^df, ( 1 ) 


with interf er omet er visibili ti es, Vh (f). given by 
(Ivan GittertI 11934 iZernikd Il938t iThompson et al.l 
12001 1 


Vbif) 


Ij A{s, f) I{s, f) Wff) dH 



sky 


A{s,f)I{sJ) 
\/l — P — m? 


Wi{f) dm, 


( 2 ) 

(3) 


where, I{s, f) and A(s, /) are the sky brightness and an¬ 
tenna’s directional power pattern, respectively, as a func¬ 
tion of frequency (/) and direction on the sky denoted 
by the unit vector s = {l,m,n), Wi(/) denotes instru¬ 
mental bandpass weights, W{f) is a spectral weighting 
function that controls the transfer function in the delay 
transform, dll = (1 — P — rnf')~^l'^ dl dm is the solid an¬ 
gle element to which s is the unit normal vector, and c 
is the speed of light, t = b ■ s/c is the geometric delay 
between antenna pairs measured relative to the zenith 

and provides a mapping to position on th e sky. _ 

The delay power spectrum is dehned as (jParsons et al.l 
Paper I): 


Pd(fcL,A:||) = |H(t)P 


\X^AB 


fD^AD\ 
[~AB~) 


with 


k^ = 


2^(1) 

D 


2-kt E{z) 
c{l + 



(4) 

(5) 

( 6 ) 


where, is the effective area of the antenna, AB is the 
bandwidth, A is the wavelength of the band center, fee is 
the Boltzmann constant, /21 is the rest frequency of the 
21 cm radiation of H i, z is the redshift, D = D{z) is the 
transverse comoving distance, AD is the comoving depth 
along the line of sight, and h, Hq and E{z) = [Hm(1 + 
z)^ -I- Hk(l + zy + are standard cosmology terms. 

In this paper, we use Hm = 0.27, Ha = 0.73, Hk = 
1 — Hm — Ha, and Elo = lOOkms”^ Mpc“^. Pd{kj_,k\\) 
is in units of K^(Mpc//i)^. 

The defining characteristics of the pitchfork signature 
are understood as follows. The steep rise in subtended 
solid angle near the horizon for a fixed delay bin size 
significantly enhances the integrated emission near the 
horizon delay limits in wide-field measurements. This is 
found to be true for diffuse emission even on wide an¬ 
tenna spacings because their foreshortening toward the 
horizon makes them sensitive to large angular scales that 
match the inverse of their foreshortened lengths. In the 
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following sections, we present an observational confirma¬ 
tion of this effect. 

3. MWA OBSERVATIONS 

The MWA instrument configuration, EoR observing 
strategy, and analysis procedure applied to individual 
snapshots used in this study are already described in 
Paper I and the references therein. In that paper, we 
analyzed two observations - an off-zenith pointing that 
included significant Galactic plane contributions and a 
zenith pointing in which the Galactic plane was signif¬ 
icantly absent. The off-zenith snapshot was useful for 
demonstrating the mapping between delay spectra and 
sky locations, establishing primary causes of foreground 
contamination, and devising a technique for mitigating 
foreground contamination. Pitchfork signatures that re¬ 
sult from low-level ubiquitous diffuse emission are fainter. 
Since Galactic plane emission, even from directions far 
away from the primary field of view, can potentially 
swamp the fainter pitchfork signatures, zenith pointings 
that have a maximum avoidance of the Galactic plane 
are preferred for this study. 

To reduce thermal fluctuations while maintaining co¬ 
herence, it is essential to average independent data sets 
obtained over the same region of sky with identical in¬ 
strument settings. Hence, we select a subset of MWA 
snapshots pointed at zenith, each with a duration of 
112 s, obtained over different nights which are aligned 
to within 72 s of each other in local sidereal time (LST) 
around a mean LST of 0.04 hr. The database contains 
14 snapshots satisfying these criteria. Two of these snap¬ 
shots, which contained amplitude and phase artifacts for 
a significant duration across different baselines, are ex¬ 
cluded from our analysis. The results of this coherent 
averaging are discussed below. 

4. RESULTS 

Figure [T] shows the delay spectra obtained from a single 
snapshot of MWA data (top), averaging LST aligned de¬ 
lay spectra from 12 individual snapshots from MWA ob¬ 
servations on different nights (middle) and from modeling 
with no thermal noise fluctuations shown for reference 
(bottom). In all panels, the foreground wedge bounded 
by horizon limits (white dotted lines) is prominent. The 
bright horizontal branch of power at t ~ 0 corresponds 
to foreground emission from the main lobe of the antenna 
power pattern pointed at the zenith. 

In the single snapshot (top), similar to the one used 
in Paper I, faint features associated with the pitchfork 
signature are visible near the horizon limits. But the 
high level of thermal fluctuations makes their significance 
marginal. In contrast, the dynamic range in the averaged 
data (middle) is a factor > 10 higher (in the delay power 
spectrum) relative to that in a single snapshot and is 
consistent with the improvement expected from averag¬ 
ing 12 independent snapshots. Hence, the foreground 
power near the horizon limits appears > 10 times more 
prominent. Also, faint horizontal features that are not 
seen in the single snapshot appear at r ~ ±0.78 /iS thus 
confirming the improvement in sensitivity. We identify 
these faint features as the response of the MWA coarse 
band edges flagged periodically every 1.28 MHz. 

In these observations, the Galactic center is just about 
to set in the west. Its signature on eastward baselines is 


10 ^ 10 ^ 10 ^ 10 '^ 10 ^ 10 ® 10 ^ ( Mpc / h )® 


[h Mpc ^ ] 

0.031 0.063 0.094 0.126 0.157 
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Fig. 1.— Delay power spectra obtained from a single snapshot 
(top), by averaging 12 snapshots of LST aligned MWA data [mid¬ 
dle)^ and from modeling with no thermal noise added (bottom). 
The fc-axis, denoted by \b\ (and kj_), represents angular (and spa¬ 
tial) scales in the plane of the sky while the p-axis, shown in r and 
/c||, denotes the spatial scales along the line of sight. White dot¬ 
ted lines are the horizon delay limits. Dynamic range in the delay 
power spectra of MWA data has increased by a factor 10 after av¬ 
eraging (middle) relative to that in a single snapshot (top). Power 
near the horizon limits caused by wide-field effects is prominent. 
Faint horizontal features at r 2 ::; ±0.78 ps are visible due to the 
effective lowering of thermal fluctuations and are the response to 
periodic coarse band edge flagging of MWA data every 1.28 MHz. 


seen in the modeled delay spectra (bottom panel) as a 
marginal brightening of the arm near the negative hori¬ 
zon limit for |b| < 125 m, consistent with our findings in 
Paper I. This spills over into higher delay modes, result¬ 
ing in the faint (< 10^ (Mpc//i)^) vertical stripes at 
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|6| < 50 m. The corresponding vertical feature is identi¬ 
fied in the averaged data as well. 

In order to show that low-level ubiquitous diffuse emis¬ 
sion is a significant contributor to the pitchfork signa¬ 
ture, contribution from any strong emission from near 
the Galactic center needs to be minimized. This is best 
illustrated with northward antenna spacings that map 
any residual emission from this region to r ~ 0 and thus 
reduce the impact on higher delay modes (Paper I). Fig¬ 
ure [2] shows the averaged delay power spectra on three 
selected baseline vectors oriented northward. Data and 
noiseless models are shown in black and red, respec¬ 
tively. The horizontal dotted black line denotes the rms 
of thermal fluctuations estimated from data. The ver¬ 
tical dashed line denotes horizon delay limits, and the 
vertical dotted-dashed lines denote delays at which the 
responses to coarse band edge flagging are expected. 


fc|| [h Mpc ^ ] 



Fig. 2.— Delay power spectra on three antenna spacings ori¬ 
ented northward, obtained by the coherent averaging of 12 snap¬ 
shots aligned in LST. The averaged data and models are shown in 
black and red, respectively. The antenna spacings are 83.5 m {top), 
104.2 m (middle), and 171 m (bottom). The horizontal dotted line 
is the rms of thermal fluctuations. The vertical dashed lines de¬ 
note the horizon delay limits. The vertical dotted—dashed lines at 
T = ±0.78 /rs correspond to the grating responses of the periodic 
flagging of bandpass at intervals of 1.28 MHz. The peaks close 
to the horizon delay limits are distinctly visible at ~ 10-1000 a 
levels. Differences between the model and the data are primarily 
attributed to uncertainties in the foreground model and the MWA 
tile power pattern. 

We focus on the prominent peaks in data near the hori¬ 
zon limits. Typically, the power near the negative hori¬ 
zon limit is seen with a signal-to-noise ratio ~ 10-100, 


while that around the positive horizon limit is ~ 100- 
1000. 

There is a remarkable agreement in broad morphology 
between the data and the model. However, some differ¬ 
ences in the amplitude scales are noted. For instance, 
the emission near the positive horizon limit is higher 
in the data than in the model in both figures. We at¬ 
tribute such differences to uncertainties in the foreground 
model, the MWA tile power pattern, thermal fluctua¬ 
tions, and other uncertainties noted in Paper 1 which 
limit a more thorough quantitative comparison between 
the model and the data. In fact, this lays further empha¬ 
sis on the need for the following - extending the footprint 
of surveys matching the frequency and angular resolution 
of observat ions such as the MWA Com missioning Survey 
(MWACS; iHurlev- Walker et al.l 1201411 and the Galactic 
and E xtragalactic MWA Survey (GLEAM; iWavth et al.l 
l2015fl to cover the entire hemisphere, and a detailed char¬ 
acterization of the power pattern over the entire hemi- 
sp here comprehensive ly extending studies such as those 
in iNeben et al.l (|2015ll . 

We note that reducing uncertainties will only change 
the relative strength of the pitchfork signature in our 
model. However, the effects giving rise to this signa¬ 
ture are generic to all wide-held measurements. Thus, 
the extremely high signihcance detection of foreground 
emission near the horizon limits is a robust conhrmation 
of the predicted effects of wide-held measurements. 

5. IMPACT ON INSTRUMENT DESIGN 

The delay spectrum maps the geometric delays to po¬ 
sitions of foreground objects on the sky. Thus the direc¬ 
tional power pattern of the antenna has a direct impact 
on the delay spectrum. In fact, from equations [U [21 and 
m the delay power spectrum scales as the square of the 
directional power pattern of the antenna. Since the con¬ 
tamination in the EoR window is strongly dependent on 
sources of emission close to the horizon, the design of 
antenna apertures that suppress their sensitivity toward 
the horizon is important. 

We approached this issue in Paper 1 via simulations of 
different antenna apertures - a dipole, a 4 x 4 phased 
array, and a dish. Among these, a dish provides the best 
quality, as evaluated from a foreground contamination 
viewpoint. With new instruments such as HERA and 
SKA on the horizon, we extend this discussion to further 
emphasize the need to comprehensively study the effect 
of their antenna apertures on foreground contamination. 

HERA will deploy 14 m dishes, with 331 of them in 
a closely packed hexagonal array and another 21 as out¬ 
liers. With the advantage of enormous redundancy and 
unprecedented sensitivity, it will address the following 
key questions - what objects first lit up the Universe and 
reionized the neutral IGM? Over what redshift range did 
this occur? And how did the process proceed, leading to 
the large-scale galaxy structure seen today? The fixed 
dishes will observe the sky drifting overhead. Based on 
simulations of uniformly illuminated circular disks, the 
response near the horizon is < —40 dB, which is over 
20 dB lower than that of the MWA. Tapering the re¬ 
flectivity of the dish is being studied to further reduce 
this response. The responses of the actual dishes will be 
tested shortly. 

With the SKA at low frequencies, deploying “stations” 
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of ^ 35 m diameter that act as aperture arrays consist¬ 
ing of 256 pseudo-randomly placed vertical log-periodic 
dipole antennas is being considered E3E3 The voltage 
beam of the station will be obtained by a phased ad¬ 
dition of the dipole responses. Since the power pattern 
of a baseline consisting of two stations will be effectively 
be a product of two voltage beams, it will homogenize the 
fluctuations in the power patterns of individual stations. 
Even with simple models of such phased array stations, 
we estimate the typical horizon response of these beams 
to be —30 dB to —40 dB compared to the zenith, yielding 
a response that is at least 10 dB better than the MWA. 
A targeted optimization of the antenna layouts of sta¬ 
tions is under active study and could yield even better 
responses. 

From these simple antenna aperture models, HERA 
and SKA should mitigate leakage from the pitchfork into 
the EoR window by ~ 40 dB relative to the MWA when 
expressed in units of power spectrum. This will signifi¬ 
cantly increase the likelihood of detecting the cosmolog¬ 
ical H I power spectrum in sensitive three-dimensional 
fc-modes adjacent to the pitchfork. The precise beam re¬ 
sponses and the resulting foreground leakage, especially 
from near the horizon, will prove to be critical and timely 
inputs to the actual aperture designs of HERA dishes 
and SKA stations, which are currently nearing their fi¬ 
nal stages. 

6. SUMMARY 

Using deeper MWA data, we have confirmed with a 
high significance the earlier prediction of a characteris¬ 
tic pitchfork morphology where wide-held EoR measure¬ 
ments suffer from signihcant foreground contamination 
from near the horizon. This has important implications 
for the instrument design and data analysis of future in¬ 
struments such as HERA and SKA. 

Careful aperture designs that suppress the response 
near the horizon, and hence also the leakage from the 
pitchfork, will signihcantly avoid contamination in k- 
modes considered critical for cosmological signal detec¬ 
tion. Precise modeling is thus required to gain a complete 
understanding of the characteristics of the cosmological 
signal and the foregrounds. 

Foreground and instrument models serve as inputs 
for power spec tr um estimation te chniques (see, e.g., 
iLiu et alll20i4al lbl: lDillon et al.ll2015ll . Input models that 
ignore wide-held effects will provide over-valued weights 
in and around the pitchfork modes, and hence result in 
sub-optimal results. Thus, conhrmation of the pitch- 
fork effect has important implications for optimum power 
spectrum estimation. 

In Paper I, we proposed a selective bagging of data 


on different baselines that can potentially mitigate fore¬ 
ground contamination by two orders of magnitude. Fol¬ 
lowing the conhrmation presented here, efforts are under¬ 
way to incorporate this proposed foreground mitigation 
technique into the MWA data analysis. 

For future work, we plan to extend our analysis to 
HERA. Based on Paper I, a dish will have a much desir¬ 
able Fourier response from a foreground contamination 
viewpoint. One of our objectives is to forecast the per- 
baseline foreground contamination as a function of LST 
in order to tune the HERA observing strategy and data 
analysis to maximize sensitivity to the EoR signal. 
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